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ABSTRACT 

We present the results from a search for HI emission from a sample of newly discovered 
dwarf galaxies in the M8 1 group. HI is detected in three galaxies, all of which are classified 
as BCDs. The HI masses of these galaxies are ^ 10^ M©, making these some of the lowest 
mass BCDs known. For these three galaxies FUV images (from GALEX) and Ha images 
(from the Russian 6m BTA telescope) are available. The Ha emission is very faint, and, in 
principle could be produced by a single O star Further, in all cases we find offsets between 
the peak of the FUV emission and that of the Ha emission. Offsets between the most recent 
sites of star formation (i.e. those traced by Ha) and the older sites (i.e. those traced by FUV) 
would be natural if the star formation is stochastic. In spite of the expectation that the ef- 
fects of mechanical feedback from star formation would be most directly seen in the smallest 
galaxies with low gravitational potentials, we only see tentative evidence of outflowing HI gas 
associated with the star forming region in one of the galaxies. 

Key words: galaxies: dwarf - radio lines: galaxies - ultraviolet: galaxies - ISM: jets and 
outflows 



1 INTRODUCTION 

Accurate determination of the faint end slope of the Luminosity 
Function in groups of galaxies is very important to resolve the dis- 
crepancy between the number of dwarf galaxies observed and that 
predicted by dark matter simulations of hierarchical structure for- 
mation l |Klypin et al. | 19991 IMoore et al. 11999 1, especially in low- 
density environments (Barkhouse et al. 2007, Roberts et al.|2007] l. 
An excellent candidate low-density environment group for carrying 
out a deep census of dwarfs is the M8 1 group of galaxies, for which 
it is possible to detect galaxies as faint as M_h ~ —9, and for which 
accurate TRGB distances can be determined for individual galax- 
ies. jCMboucasT^^achents^T present a census of 
the dwarfs in this group, determined from a large area survey us- 
ing the Megacam instrument on the CFHT. A total of 22 candidate 
dwarf galaxies were found in this survey. Here we present follow up 
HI observations of some of these newly discovered dwarf galaxies. 

The observations were done at the Giant Metrewave Radio 
Telescope (GMRT) and targeted a sub-sample of galaxies that 
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showed evidence for active star formation, and as such was ex- 
pected to have a significant HI content. At the redshift of the M81 
group, HI emission from group members is expected to overlap in 
velocity with emission from HI in our Galaxy. An interferometric 
survey is hence the only reliable method for robustly detecting gas 
associated with the dwarfs in this group. Deeper follow-up observa- 
tions were obtained using the GMRT and the Westerbork Synthesis 
Radio telescope (WSRT) for the galaxies in which HI was detected. 

The target galaxies were selected based on the original Mega- 
cam CFHT survey data. However, subsequently FUV data from the 
GALEX satellite as well as Ha data from the Russian 6m BTA 
telescope also became available. We present and analyse this data 
in conjunction with the HI data from the GMRT and the WSRT. The 
galaxies for which we do have GMRT detections have HI masses 
~ 10® M0, and star formation rates ~ 10^'* Mq/jv. They repre- 
sent some of the smallest mass star forming galaxies known. Feed- 
back into the ISM from star formation is expected to be most pro- 
nounced for small galaxies, and we use our data to explore this 
issue too. 

Below, in Sec. [2] we present the sample, briefly describe the 
observations and results. The results are discussed in detail in 
Sec.|3] and the summary and conclusions presented in Sec.|4] 
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Table 1. The sample 



Name 


Q (J2000) 


5(J2000) M,/" 


axb'' Dho'= 


Type"* 




Sub/ BTAfl 
rad rad 




(h m s) 




(' X ') (') 




(Mpc) 


(kms-1) (kms-i) 



d0926+70 


09 26 27.9 


+70 30 24 


-9.7 


0.8x0.6 


0.40 


dl 


T T+0.2 

-'■-'-0.2 






d0944+71 


09 44 34.4 


+71 28 57 


-12.5 


1.3x0.8 


0.90 


dl/dSph 


3 4+0.1 






d0958+66 


09 58 48.5 


+66 50 59 


-13.2 


1.6x0.8 


0.88 


BCD 


-, O+O.l 

^■^-0.1 


+33±94 


+90±50 


dl009+70 


10 09 34.9 


+70 32 55 




0.6x0.5 


0.73 


dl 








dl012+64 


10 12 48.4 


+64 06 27 


-13.4 


1.4x0.9 


0.90 


BCD 


T 7+0.1 
^■'-0.1 




+150±50 


dl014+68 


10 14 55.8 


+68 45 27 


-9.4 


0.4x0.4 


0.20 


dSph 


n n+0.3 

^■»-0.3 






dl028+70 


10 28 39.7 


+70 14 01 


-12.4 


1.4x0.9 


0.76 


BCD 


T O+O.l 

^■°-a.i 


-90±79 


-114±50 


dl041+70 


10 41 16.8 


+70 09 03 


-9.3 


0.8x0.4 


0.48 


dl 


-, 7+0.2 
-^■'-0.3 







" estimated extinction corrected r' absolute magnitude, from |Chiboucas, Karachentsev, & TullyH2009) corrected for the measured distance 
^ 28 mag/arcsec^ diamaters of the Megacam images 
Holmberg diameter 

as listed in |Chiboucas, Karachentsev, & TullyH2009t 

extinction conected TRGB distances using HST, from |Chiboucas et al.|(20T2) 
/ f Optically measured radial velocities from the Subaru telescope or 6m BTA telescope 
* Not in M8 1 group, no distance measured. 



Table 2. Observation Details 



Name 


Array 


Date of Observation 


Flux calibrator 


Phase calibrator 


Time on source 
(hours) 


Channel width 
(kms-i) 


Synthesised beam 
(arcsec^) 


Noise 
(mJy) 


d0926+70 


GMRT 


04 Sep 07, 07 Sep 07 


3C147,3C286 


0834+555.1035+564 


~6.9 


3.3 


33x32 


2.0 


d0944+71 


GMRT 


29 Oct 07 


3C147,3C286 


0834+555 


~6.1 


3.3 


38x28 


1.9 


d0958+66 


GMRT 


04 Sep 07, 07 Sep 07 


3C147,3C286 


0834+555,1035+564 


~6.0 


3.3 


34x33 


1.4 




GMRT 


01 Jan 12 


3C48,3C286 


0834+555,1035+564 


~9.2 


3.3 


40x40, 12x11 


1.6,0.9 


d 1009+70 


GMRT 


04 Sep 07, 07 Sep 07 


3C147,3C286 


0834+555.1035+564 


~6.1 


3.3 


35x32 


2.2 


dl012+64 


GMRT 


28 Oct 07, 17 Jan 08 


3C147,3C286 


0834+555 


~6.5 


3.3 


33x30, 13x13 


2.0, l.I 




WSRT 


16 Sep 11 


3C48,3C286 


3C48,3C286 


~11.9 


3.1 


42x42,23x23 


1.2,0.9 


dl014+68 


GMRT 


26 Nov 07 


3C147,3C286 


0834+555 


~7.1 


3.3 


48x43 


2.4 


dl028+70 


GMRT 


28 Aug 07 


3C286 


1035+564 


~3.1 


3.3 


39x29 


1.3 




GMRT 


02 Jan 12 


3C48,3C286 


1035+564.1313+675 


~8.0 


3.3 


43x43 , 14x14 


1.7 , 1.0 


d 104 1+70 


GMRT 


28 Aug 07, 08 Sep 2007 


3C147,3C286 


0834+555 


~9.3 


3.3 


22x18 


0.9 



Table 3. Results from HI observations 



2 SAMPLE, OBSERVATIONS AND RESULTS 



Name 


Velocity" 


Velocity'' 






(kms-i) 


(kms-i) 


(IO6M0) 


d0926+70 


-208 to 567 




< 0.31 


d0944+71 


-137 to 637 




< 0.31 


d0958+66 




89.6±0.9 , 29±2 


0.73±0.O8* 


dl009+70 


-208 to 567 






dl012+64 




189.7±0.4, 19±1 


1.04±0.09* 


dl014+68 


-139 to 636 




<0.48 


dl028+70 




-68.9±0.5 , 19±1 


1.3±0.1* 


dl04I+70 


-208 to 567 




< 0.17 



" velocity range searched in case of non detections 
^ Vpeaki ^Vso in the case of detections 

3cr upper limits in the case of non detections 
* error on mass includes flux and distance measurement errors. 



The sub sample chosen to be observed with the GMRT is listed 
in Table [T| The main selection criterion was that the galaxy should 
have been classified as a dirr or BCD, so that it was likely to contain 
HI. At the time of preparing the observing sample, dl014+68 was 
suspected to be a transition galaxy and was hence also added to the 
sample. Further, subsequent to the GMRT observations, dl009+70 
was found not to belong to the M81 group and hence no distance 
for it was measured. 

The details of the initial GMRT observations, as well as the 
later GMRT and WSRT observations of three galaxies in the sam- 
ple, are given in Table[2] Please note that the noise value quoted in 
Column (9) is the rms noise in a line free channel of the (continuum 
subtracted) data cube made at the corresponding spatial resolution 
mentioned in Column (8). The data were reduced using the As- 
tronomical Image Processing System (AIPS) using standard spec- 
tral line data analysis procedures. Total intensity and velocity field 
maps were created using the AIPS task MOMNT, with the param- 
eters in the adverb FUNCTYPE set to 'HG' and the adverb CELL- 
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Table 4. Derived parameters for detected galaxies 



Galaxy 


Dhi 
(') 


IHI 

n 


Mdyn 
(108 Mq) 


Dm 

Dho 




Mhi 


SFRpuv 
(Mq yr-1) 


THI 

(Gyr) 


(ergs s-^) 


SFRhq 

(Mq yr- 


1) 


d0958+66 


1.5 


~29 


~1.69 


1.7 


-11.90 


0.08 


5.0x10-4 


1.5 


3.7x10^6 


3.0x10- 


5 


dl012+64 


2.0 


~25 


~1.26 


2.2 


-12.24 


0.09 


6.2x10-'' 


1.8 


2.45 xlO^'^ 


1.94x10" 


-4 


d 1028+70 


2.3 


~44 


~0.54 


3.0 


-11.70 


0.17 


7.0x10""' 


1.8 


3.03x10^^ 


2.40x10" 


-4 




50 100 150 

HeliocGntric velocity (km s'^) 



0.05 
0.045 



0.02 
0.015 

0.01 
0.005 


-0.005 
-0.01 



■ d1012-i-64 










\ 






\ 
\ 






V ■ 










1 50 200 
HGliocentric velocity (km s"^} 



-100 -50 
Heliocentric velocity (km s"') 



Figure 1. HI spectra of the three detected BCDs derived using the follow-up observations, from the low resolution (~ 40 ) data cubes. The single Gaussian 
fit to each spectrum is also shown as dashed lines. 
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Figure 2. The upper row shows the overlays of the low spatial resolution HI images (contours) of the BCDs d0958+66 [A], dl012+70 [B], and dl028+64 [C] 
on their respective R band Megacam images (greyscales). The resolutions are 40 , 42 and 48 for d0958+66, dl012+70 and dl028+64 respectively (beams 
also shown as circles in the images). The contours start at the 5cr sensitivity levels for each galaxy, with each consecutive contour \/2 times the previous 
one (please see Section|2]for details about how moment maps were constructed). The lower row shows the integrated velocity fields for the three galaxies 
(d0958+66 [D], dl012+70 [E] and dl028+64 [F]) at the same spatial resolution as the images in the upper row. For each galaxy, the images in the upper and 
lower rows are of the same spatial extent. 
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SIZE set to 3 6 (Except for the highest resolution moment maps of 
dl012+64 for which the parameters were set to 3 3). 

Table |3] gives a summary of the basic deductions one can 
make from the HI observations of these galaxies. Mass limits based 
on the measured noise are provided for the non-detections. These 
are 3a mass limits assuming the galaxy to be unresolved in a 
~ 40 beam, and having a typical velocity width of 20 km s~^ tee-' 
|gum et aT. 2008 1. HI was detected in three of the sample galaxies, 
and deeper observations of these three galaxies were later done with 
the GMRT and the WSRT. We use these more sensitive observa- 
tions for all subsequent analysis. The peak and velocity spread of 
the detected HI line is given for these three galaxies, and the values 
are based on the low resolution data cubes. For each galaxy, only 
channel fluxes more than two times the rms of the fluxes measured 
in line free channels were considered while obtaining the total flux. 
As can be seen from the table, HI emission has been detected only 
from those galaxies that are classified as BCDs. The HI masses of 
the detections are all ~ IO^Mq, making these the lowest HI mass 
BCDs with interferometric observations. The integrated spectra for 
the detected galaxies are shown in Figure [T| 

Figure [2] shows an overlay of the low resolution (~ 40 ) 
GMRT and WSRT HI images of the three detected BCDs and 
the Megacam R band images. The moment and moment 1 maps 
shown for all the galaxies were constructed using the previously ex- 
plained method after applying a flux threshold of 2 times the noise 
in the line free channels. It is apparent that the HI discs extend 
much more than the stellar disc (see Column (5) in Table [4]l, and 
the stellar disc is approximately coincident with the central den- 
sity peak of HI for all the galaxies. Looking at the velocity fields 
for the galaxies, we see that in all cases there is an overall velocity 
gradient approximately aligned with the morphological major axis. 
However, the isovelocity contours are not regular, and the velocity 
fields cannot be modelled as simple rotating disks. Not surprisingly, 
attempts to fit tilted ring models to the velocity fields did not lead 
to meaningful fit parameters. 

Table |4] lists various measured global parameters of the three 
detected BCDs. The entries in the table are: Column(l) galaxy 
name, Column(2) HI diameter Dhi, from the 10^^ cm~'^ column 
density contour measured in the low resolution Moment map 
(shown in Figure |2|. For all the three galaxies, this contour was 
near or above the 3(7 sensitivity contour. Column(3) an indica- 
tive inclination angle (i) of the HI disk. The inclination angle was 
determined by measuring the axial ratio of the above mentioned 
10^® cm~'^ column density isophote and assuming an intrinsic ax- 
ial ratio of 0.57 ( [Roychowdhury et al.|2010^ if one assumes that 
the HI discs are oblate spheroids. The inclination angle estimated 
in such a manner is highly uncertain, and is only used here to quan- 
tify what is apparent from looking at the HI discs of the galaxies: 
d0958-l-66 and dl012-l-64 have very face-on HI discs and dl028-l-70 
has a slightly inclined disc. Column(4) an indicative dynamical 
mass estimated as Mdyn = (Al/50/2 sin(i))^ x {Dm/2)/G. Col- 
umn(5) the ratio of the HI to the Holmberg diameter. Column(6) 
is the absolute blue magnitude obtained by applying the distance 
measurements of the galaxies to the blue magnitudes estimated by 
eye from POSS-II images tabulated in Chiboucas, Karachentsev,, 
|& Tully| (2009 l. Column(7) the HI mass-to-light ratios of the three 
galaxies, using the magnitudes tabulated in column (6). Column(8) 
gives the total SFR as deduced from FUV emission (SRFpuv) as 
measured by GALEX, and using the calibration given in Kennicutt| 
( |1998j l (more details on this are given below). Column(9) The gas 
consumption timescale thi defined as thi = Mhi/SFR^'^^. Col- 



umn(lO) The total Ha luminosity as determined from the 6m BTA 
telescope. Column(l 1) the SFR computed from the Ha luminosity, 
using the calibration in |Kennicutt| ( [ 1 998) (more details on the last 
column are also given below). 

[Gil de Paz, Madore & Pevunova| ( [2003) suggest a set of quan- 
titative criteria for classifying a galaxy as a BCD. Unfortunately 
these involve photometry in the B, R and K bands, which is not 
available for our galaxies. However from the measured Ha lu- 
minosity and Mb for our galaxies we see that they follow the 
same general trend as the brighter BCDs in [Gil de Paz, Madore| 
|& PevunovaH2003| (see their Fig. 4e). Also, the gas consumption 
timescales for our galaxies is similar to that typical for brighter 
BCDs ( [Sanchez Almeida et al.|2008) . It is however interesting to 
note that the relative content of neutral hydrogen in these galax- 
ies is lower than that seen for brighter BCDs with HI observations, 
more so keeping their faint blue magnitudes in mind (e.g., compar- 
ing with Figure 2 in [Huchtmeier, Krishna & Petrosian|2005| the 
ratios are all below the best fit to their data). 

The SFR from the FUV was computed by taking the to- 
tal FUV flux contained within the best fit ellipse to the 28th 
magnitude/arcsec^ isophote of the CFHT Megacam images. The 
centre and principal axes of this ellipse are given in Table [T] This 
FUV flux was then corrected for the local extinction expected from 
dust in our galaxy. Note that the presence of foreground galactic 
cirrus with substantial fine-scale structure in the line of sight to- 
wards the M81 group makes this correction somewhat uncertain 
( [Sandage|1976[ [Pavies et al.[[2010) . No correction was made for 
extinction internal to the BCDs, since this is expected to be small. 
The corrected FUV flux was converted to a SFR using the stan- 
dard calibrations for obtaining SFRs given in Kennicutt[ ( [T998l l. 
Note that this calibration assumes solar metallicity and a constant 
star formation rate over the last ~ 100 million years. From the 
luminosity-metallicity relation one would expect that our sample 
galaxies would have a lower than solar metallicity (see below), and 
further that since they are BCDs they may be undergoing a recent 
starburst. However, in view of the fact that there is some uncer- 
tainty in the foreground FUV extinction, no correction was made 
for these two effects. Similarly, the SFR was computed from the 
Ha flux assuming the standard calibration from |Kennicutt[ ( |T99 8| l . 
As discussed further below, this is unlikely to be a reliable estimate 
for the SFR in these galaxies. 



3 DISCUSSION 

3.1 HI non-detections 

The newly-discovered dwarf galaxies in the M81 group ( [Chibou-j 
[cas, Karachentsev, & Tully|2009[ l chosen for follow-up HI obser- 
vations were the brightest ones, and were thought to be star forming 
galaxies as complexes of young blue stars seemed to be present in 
all of them on initial visual inspection. Subsequent Hubble Space 
Telescope (HST) observations revealed that d 1009-1-70 is a back- 
ground galaxy not in group (possible companion to UGC 5423 
at a distance of 8.7 Mpc) and hence its HI emission may be too 
faint to be detected. Of the other five galaxies with no detected HI 
emission, follow-up photometry using HST data has revealed that 
d0944-l-71 (initially could not unambiguously be classified as ei- 
ther dl or dSph), dl014-l-68 (initially thought to be transition type 
dwarf) and d 104 1-1-70 (initially classified as dl as mentioned in Ta- 
ble[TJ are all dSph galaxies. Even if the above galaxies have young 
stellar complexes, these are presumably small (as can be said from 
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Table 5. Revised SFRs from IGIMF theory 



Galaxy 


(M0 yr-l) 


(M0 yr-l) 


SFR?'<1 

Ha 

(M© yr-l) 


SFRSi"i 
(Mq yr-l) 


d0958+66 


1.7x10-3 


9.1x10-1 


1.7x10-3 


9.9x10-* 


dl012+64 


2.0x10-3 


1.1x10-3 


4.3x10-3 


2.3x10-3 


dl 028+70 


2.2x10-3 


1.2x10-3 


4.8x10-3 


2.5x10-3 



+66°51'20' 




+50'.)0" 



+50'30"| 



9''5a'"54' 52' 50' 48' 



10"12"'50' -18' 
Right Ascension (J2000) 



40" 35' 



Figure 3. Comparison of HI, FUV, Ha, and optical data for the BCDs d0958+66 [1], dl012+70 [2], and dl028+64 [3]. The first row [a] shows the highest 
spatial resolution HI contours in cyan (with beam sizes 12 X 11 , 13 X 13 and 14 X 14 for d0958+66, dl012+70 and dl028+64 respectively, which 
are also shown as circles) overlayed on the FUV emission smoothed to 6 resolution shown in greyscale. The contours start at the 3a- sensitivity levels for 
each galaxy (please see Sectionplfor details about how moment maps were constructed) and subsequent contours are at intervals of \/2. The second row [b] 
shows Ha emission contours in red overlayed with FUV emission contours in blues (both smoothed to 6 resolution, beam sizes also shown as circles). For 
d0958+66, the contours start at the 3(j sensitivity levels for FUV and Ha emission, and are spaced apart by \/2. For dl012+70 and dl028+64, the contours 
start at the 5cr sensitivity levels for FUV and Ha emission, and are spaced apart by 2. The greyscale images of FUV in row [a] start at the level of the first blue 
contours in row [b] . The lower row [c] shows the RGB images constructed from HST observations, for the same spatial region as the images in the top row, 
separately for each galaxy. In the RGB images, the red and green channels are the 1 (F814W) and V (F606W) band images, whereas the blue channel image is 
constructed as 2x V-I. For d0958+66 and dl028+70 the images are from WFPC2, whereas for dl012+64 the images are from ACS. 
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Figure 4. The HI emission in different velocity cliannels of dl012+64, imaged at a spatial resolution of 23 (shown as circles), is shown as blue contours. 
All cutouts represent the same spatial area with the central heliocentric velocity for each channel mentioned in the top right comer (in km s^^). The HI (blue) 
contours start at the 3(T sensitivity level, and are spaced apart by 1.2. The HST ACS V (F606W) band image is shown as greyscale in all images, and represents 
the stellar disc. The grey contours shown in the last image is the integrated velocity field of the galaxy at the same resolution at which the HI emission for 
each channel is shown (please see Section[2]for details about how moment maps were constructed). The velocity which corresponds to any particular contour 
is mentioned (in km s~^). 



their faintness) and are being fed by HI reservoirs with masses be- 
low the detection limit of the present observation. Follow-up HST 
photometry also shows that d0926-l-70 is indeed a dl galaxy, though 
with very little star formation which explains the non-detection of 
HI using similar reasoning as stated above for the dSph galaxies. 
For all the non-detections, there is the possibility that they have 
systemic velocities outside the velocity range searched. But this 
is an unlikely possibility considering the fact that TRGB distance 
measurements constrain them to be in the M8 1 group, and the mea- 
sured velocities of other M81 group members jKarachentsev et aT] 
|2002[ > are in the range of velocities searched. 

3.2 Hq and FUV emission 

In Figure |3] the Ha and FUV emission are shown with respect to 
the corresponding high resolution HI maps for the three BCDs. The 
moment maps shown for all the galaxies were constructed us- 
ing the previously explained method after applying a flux thresh- 
old of 2 times the noise in the line free channels for d095 8-1-66 
and d 1028-1-70, and 1.5 times the noise in the line free channels 
for d 10 12-1-64 (for which old GMRT data was used to make the 



high resolution map as WSRT cannot map at this resolution). Also 
shown are the RGB images of the three galaxies constructed from 
HST dat^ A few things are apparent from the figures: (i) although 
the star forming disc is at the central peak of the HI profile when the 
HI is imaged at a low resolution (Figure[2]l, the peaks of FUV and 
Ha emission don't concur with the peaks of the HI emission when 
HI is images at high resolution, (ii) There are offsets between the 
peaks of the Ha emission and that of the FUV emission, these are 
most pronounced in d0958-(-66 and dl028-l-70. (iii) The stellar disc 
as appears in the RGB images is better traced by the FUV emission. 

Before examining the differences between the Ha and the 
FUV emission it is useful to take a look at the total integrated fluxes 
from the entire galaxy (Table|4j. As can be seen the total Ha fluxes 
are very low - in fact they are consistent with that expected from 
a single O star. The ionized flux of O stars is critically dependent 

^ The HST images and other materials for the three BCDs can be found at 
http://edd.ifa.hawaii.edu Select the catalog CMDs/TRGB and then select 
the desired object in the column Name/CMD. Alternate names for the galax- 
ies are d0958+66 = PGC0028826, dl012-i-64 = PGC0029735, dl028-H70 = 
PGC5056941. 
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Figure 5. Separate spectra for four different regions in tlie galaxy dl 01 2+64. In the top left hand comer the moment map in HI is shown in grey scale, with 
the four regions marked in different colours. The spatial extent of the image and the resolution at which HI has been mapped (also shown as a circle) are 
identical to the images shown in Figure|4] The spectra for each of the four regions are shown in the same colour in which the regions have been marked. The 
dashed line in each spectrum is at 3 times the rms noise of the fluxes in the line free channels for that particular spectrum. The red spectrum is for the region 
which shows the suspected HI outflow, marked here by the grey circle around the concerned channels. 



on the metallicity. Since we do not have abundance quality optical 
spectra to determine the metallicity, the metallicity of our galaxies 
can be estimated from the L-Z relation. The absolute blue mag- 
nitudes of our BCDs are ~ —12.0 mag. From the L-Z relation 
for BCDs (see |Ekta & Chengalur|20f0| l, the expected metallicity 
is 12+log(0/H) = 7.30, i.e. about 0.05 times the solar metallicity. 
From the 0.05 solar metallicity model atmospheres given in |Smith,| 
|Norris, & Crowther|j2002^ , the Ha flux expected from an ioniza- 
tion bounded HIT region excited by a single 04 star is 3.4 x 10^^ 
ergs/s, the corresponding number for an 08 star is 3.4 x 10^® ergs/s. 
These values bracket the measured Ha fluxes for our galaxies. It is 
worthwhile to note that the L-Z relation for BCDs used here is quite 
well constrained, as it leaves out the deviant extremely metal defi- 
cient (XMD) BCDs (see their Figure 2). Even if the galaxies con- 
sidered here turn out to be XMD BCDs, then the fluxes expected 
from the HII regions will only increase. 

From the SFRs deduced using the total FUV and Ha fluxes 
listed in Table|4]one can also see that for all the three BCDs the de- 
duced SFR^"^ is > SFR^^". This is typical for low mass galaxies 



I Lee et el. | 2009[[Hunter, Elmegreen & Ludk a 20 1 ; Roychowdhury 
|et al.|201 1 1. As discussed in detail in these (and other papers), this 
is likely due to the fact that at low star formation rates, the num- 
ber of massive stars that are present at any given time is subject to 
substantial statistical fluctuations. The Ha emission (which comes 
from massive M> 17 M© stars) will in general be weaker than that 
predicted from the FUV luminosity (which comes from intermedi- 
ate M > 3 Mq stars). 

A lower Ha flux compared to what would be expected from 
the observed FUV emission and a standard Salpeter IMF is also 
predicted by the integrated galaxial initial mass function (IGIMF) 
theory ( [Weidner & Kroupa|2005^ . In this theory, the clustered na- 



ture of star formation, along with the assumption that the most mas- 
sive star formed in a cluster depends on the total mass of the cluster, 
results in the total sampled IMF to differ from the Salpeter IMF. 
The theory predicts that the FUV and Ha luminosities decrease 
non-linearly with the total SFR; the effect is most pronounced 
for galaxies with low SFRs. [Pflanun-Altenburg & KroupaH2009| l 
show that recalibrations using the IGIMF explain to a large ex- 
tent the drastic decrease in Ha luminosities for dwarf galaxies. Ta- 
ble |5] shows the recalibrated SFRs for our BCDs, from both FUV 
and Ha, for two sets of IGIMF models (i.e. those called standard 
and minimal- 1 in IPflamm-Altenburg, Weidner & Kroupa| p007| l) 
which cover the full range of IGIMF models. Colunms(2) and (3) 
are the recalibrated SFRs from FUV using calibrations given in 
[Pflamm-Altenburg, Weidner & Kroupa|p009p , while Columns(4) 
and (5) are the recalibrated SFRs from Ha using calibrations given 
in 'Pflamm-Altenburg, Weidner & Kroupa (20071. The recalibrated 
SFR estimates match quite well for d0958+66, but for the other 
two BCDs, the SFR as derived from Ha is approximately twice the 
SFR estimated using FUV flux. From the above discussion, it is not 
clear whether there is any way to distinguish between the effects of 
stochasticity in high mass star formation and variable IMF on the 
Ha and FUV SFRs. But recent modelling using spectral synthe- 
sis codes of star clusters which include stochastic sampling of the 
IMF (e.g., [Fumagalli, da Silva, & Krumholz|2011[|Eldridge|20T2l l 
suggest that in the low SFR regime it is the determining factor. 
The data presented in 'Fumagalli, da Silva, & Krumholz (2017) (see 
their Figure 1) is for starbursts with slightly higher FUV fluxes than 
in the current case, however our data is consistent with the trends 
shown in their figure. As those authors emphasize, stochastic sam- 
pling of the IMF could occasionally result in small clusters with a 
disproportionate number of massive stars (which otherwise would 
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be evidence for a "top heavy" IMF), and this appears to be the case 
here, especially for the galaxy d0958+66. 



of neutral gas outflow associated with the star forming region of 
dl012+64. 



3.3 Signs of feedback from star formation 

We expect that the effects of mechanical feedback from star for- 
mation into the ISM would be more pronounced and easier to de- 
tect in smaller low mass galaxies. For example, [Mac Low & Fer-| 
|rara| ( [1999; ) modelled the effects of supernova explosions from star- 
bursts on the interstellar medium of dwarf galaxies to show that for 
galaxies with total gas mass ^ 10^ M© a significant fraction of 
the ISM could be expelled. This result however is sensitive to both 
the concentration of the starburst and the geometry of the ISM. 
Feedback from star formation is also believed to play an important 
role in determining the final dark matter distribution. For example, 
|Govemato et al.| ( [20I0^ show that supernova driven outflow of low 
angular momentum gas from the centres of star forming galaxies 
inhibit the formation of bulges and can account for the shallow (in- 
stead of cuspy) central dark matter profiles seen for nearby gas rich 
dwarfs. This weakening of the central dark matter cusp occurs even 
if there is no escape of baryonic matter from the galaxy. On the ob- 
servational front, [van Eymeren et al.| ( [20I0| > as well as [Cannon et| 
[al. ( 201 1 1 both find evidence for feedback from star formation into 
the neutral ISM of dwarf galaxies. But even though our galaxies 
represent some of the smallest mass, actively star forming systems, 
evidence for feedback between the star formation on the neutral gas 
distribution in two of them (d0958-l-66 and dI028-l-70) is not clearly 
evident. 

The only evidence we see of mechanical feedback form star 
formation on the HI distribution is in the galaxy dIO 12-1-64. Fig- 
ure|4]shows the HI channels maps for this galaxy at an intermediate 
resolution (chosen to optimise between localising the HI emission 
to a particular region and having enough signal-to-noise ratio to 
pick up the emission). The stellar disc, and the integrated velocity 
field at the same resolution (moment 1 map) are shown for refer- 
ence. The velocity field shows a South-East to North alignment as 
one moves forward in the velocity space. The channels of interest 
are the last two ones shown in the Figure, with central heliocen- 
tric velocities of 211.5 km s~^ and 214.7 km s~^. Here the HI 
suddenly moves out of the South-East to North flow for increas- 
ing velocity, and returns to the central region of the HI disc of the 
galaxy. The noteworthy fact is that the HI emission in these two 
channels is right at the edge of the stellar disc, and overlapping 
with the Ha emission coming from recent star formation (compare 
with Figure 3, 2b and 2c). To confirm that there is high velocity 
HI concentrated in only one region, we plot the HI spectrum for 
four different regions around the peak of HI emission in Figure|5] 
There is statistically significant detection of HI in two channels at 
heliocentric velocities above 210 km s~^ in only the region to the 
North-East, as is seen from Figure [4] too. The total mass of HI in 
the relevant channels is ~ 1.4 x lO'^M©. Taking the mean ve- 
locity difference between the channel velocities and the systematic 
velocity of HI ( ~ 23.5 km s~^), the estimated kinetic energy being 
carried away by the gas is hence Eke ~ 8 x lO*'"* ergs. Consid- 
ering the wind mechanical energy output from an OS star within a 
typical lifetime of ~ 3 Myr is ~ 1.1 x 10''* ergs (Smith, Nor- 



|ris, & Crowther|2002[ l, just the mechanical energy output from the 
Ha emitting cluster in dl0I2-l-64 may not be enough to produce 
the gas outflow seen. But this much of energy can easily be ac- 
counted for even by the typical energy output of a single supernova 
(i.e. ~ 10^^ erg). Hence, what we see is most probably signature 



4 SUMMARY AND CONCLUSIONS 

Our search for HI in the newly discovered dwarf galaxies in the 
M81 group has resulted in three new detections. The detected 
galaxies are some of the lowest mass BCDs known, with star for- 
mation rates ~ 10"'' M0/yr. The star formation can be traced 
by both Ha and FUV emission, and as is the case for most other 
dwarfs, the observed Ha flux is less than what would be predicted 
from the star formation rate deduced from the FUV. This is likely 
to be a result of stochastic sampling of the upper end of the mass 
function. The Ha emission is also offset from the peak of the FUV 
emission. We find tentative evidence of HI gas outflow at the edge 
of the stellar disc of one of the galaxies. Even though this result 
agrees with the expectation that the effects of mechanical feedback 
from star formation would be most directly seen in the smallest 
galaxies, surprisingly no clear evidence of feedback associated with 
star formation is seen in the HI discs of the other two galaxies. 
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